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We consider the possibility that the heavier CP-even Higgs boson (H^) in the minimal supersym¬ 
metric standard model (MSSM) decays invisibly into neutralinos in the light of the recent discovery 
of the 126 GeV resonance at the CERN Large Hadron Collider (LHC). For this purpose we consider 
the minimal supersymmetric standard model with universal, non-universal and arbitrary boundary 
conditions on the supersymmetry breaking gaugino mass parameters at the grand unified scale. 
Typically, scenarios with universal and nonuniversal gaugino masses do not allow invisible decays 
of the lightest Higgs boson (h°), which is identified with the 126 GeV resonance, into the lightest 
neutralinos in the MSSM. With arbitrary gaugino masses at the grand unified scale such an invisible 
decay is possible. The second lightest Higgs boson can decay into various invisible final states for a 
considerable region of the MSSM parameter space with arbitrary gaugino masses as well as with the 
gaugino masses restricted by universal and nonuniversal boundary conditions at the grand unified 
scale.The possibility of the second lightest Higgs boson of the MSSM decaying into invisible channels 
is more likely for arbitrary gaugino masses at the grand unified scale. The heavier Higgs boson decay 
into lighter particles leads to the intriguing possibility that the entire Higgs boson spectrum of the 
MSSM may be visible at the LHC even if it decays invisibly, during the searches for an extended 
Higgs boson sector at the LHC. In such a scenario the nonobservation of the extended Higgs sector 
of the MSSM may carefully be used to rule out regions of the MSSM parameter space at the LHC. 

PACS numbers: 12.60.Jv, 14.80.Da, 14.80.Ly, 14.80.Nb 


I. INTRODUCTION 

With the discovery [1-4] of a neutral state around a mass of 126 GeV at the CERN Large Hadron Col¬ 
lider (LHC), a new era has begun in our quest for the understanding of the fundamental constituents of matter 
and the forces between them. Although, the properties of the discovered state are consistent with the properties 
of the standard model (SM) Higgs boson, it also opens up a window for searches of new physics. It is well 
known that the SM Higgs sector suffers from the naturalness and hierarchy problems, thereby rendering a light 
Higgs boson technically unnatural. The most popular, and well motivated, extension of the SM which renders a 
light Higgs boson technically natural include supersymmetric models [5], of which the minimal supersymmetric 
standard model (MSSM) is perhaps the most economical, and hence compelling [6]. In the MSSM the Higgs 
spectrum is richer as compared to the Higgs sector of the SM, and consists of two Higgs superfields {Hi and H 2 ). 
After spontaneous symmetry breaking the model contains two CP even Higgs bosons (ft.°, M^o < Mho), 
one CP odd Higgs boson (A°), and two charged states (H^). At the tree level the Higgs sector of MSSM is 
rather constrained and is described by two parameters, usually taken to be the mass of A° (Ma) and the ratio 
of the vacumm expectation values of the two Higgs fields, tan/3 = {H 2 )/{Hi). Discovery of more than one 
Higgs boson at the LHC would point towards an extension of the SM, of which supersymmetry as embodied 
in the minimal supersymmetric standard model is a leading candidate. It would, therefore, be important to 
discover an extended Higgs sector at the LHC, if it exists. However, from the point of view of supersymmetry, it 
is crucial to discover the supersymmetric partners of the SM states as predicted by the MSSM such as squarks, 
gluinos and sleptons, as well as neutralinos and charginos. In the absence of any signal for supersymmetric 
particles, it would, therefore, be appropriate to ask the question whether non observation of an extended Higgs 
sector would imply that the new physics is at a much higher scale. This question is intimately connected with 
the decay patterns of the Higgs bosons of the MSSM. It is possible that the Higgs bosons of the MSSM may 
decay into some of the lighter particles of the supersymmetric spectrum at a rapid rate. This, in particular, 
would include light neutralino pairs. This decay can be an important decay channel in certain regions of 
the MSSM parameter space. Hence we need to study the invisible decay of the second lightest Higgs boson 
extensively. Another aspect of this issue is that even though MSSM could be ruled out, there are appealing 
alternatives, namely the nonminimal or next-to-minimal supersymmetric standard model (NMSSM) whose 
Higgs [7-14] and neutralino [15-17] sectors are richer than the MSSM, thereby increasing the possibility of an 
invisibly decaying Higgs boson. Furthermore, in the NMSSM there is a possibility of a low mass pseudoscalar 
in the Higgs spectrum. Thus, there is a distinct possibility of the scalar Higgs boson(s) decaying invisibly into 
these pseudoscalars [18]. 

The invisible Higgs decay width has been constrained by various groups by performing fits of the signal 
strengths in various search channels using the latest LHC Higgs data [19, 20]. Direct searches for invisible 
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decaying Higgs produced in association with a Z boson has been carried out by the ATLAS [21, 22] and 
CMS [23] collaborations at the LHC and they have excluded branching ratio of more than 65% and 75%, 
respectively, with 95% CL. The CMS collaboration [24] has also carried out a similar search for invisible 
branching ratio of the Higgs boson produced in the vector boson fusion process and put an upper limit of 69% 
on the invisible branching ratio of the lightest Higgs boson. All these searches have in turn put constraints 
on the MSSM parameter space. Although specific regions of the MSSM parameter space have been ruled out 
by these experiments, there is still a large portion of the MSSM parameter space which remains unexplored. 
Recently, the question of the invisible decays of the lightest Higgs boson in the context of the MSSM and the 
NMSSM have been discussed in detail [18, 25]. In this paper we carry this investigation further by analyzing 
the decay patterns and the invisible decays of the heavier Higgs boson (i7°) of the MSSM. For this purpose 
we shall identify the state observed near 126 GeV at the CERN LHC with the lightest Higgs boson of 
the MSSM. We shall systematically study the scenarios under which such a possibility can arise and discuss 
different aspects of these scenarios. 

The plan of the paper is as follows. In Section H we discuss the relevant features of the Higgs sector of 
the MSSM, and enumerate the two regions, the decoupling and the non-decoupling regions, of its parameter 
space. In Section HI we discuss the neutralino sector and the neutralino mass matrix of the MSSM with two 
different kinds of grand unified theory (GUT) boundary conditions on its parameters, namely universal and non- 
universal boundary conditions and the composition of neutralinos. In Section IV we summarize the analytical 
expressions for the decay of second lightest Higgs boson to neutralinos in the MSSM with an appropriate 
choice of the parameter space, taking all the experimental contraints into account. In Section V we present our 
calculations and numerical results for the invisible decay of the second lightest Higgs boson and comment on 
how some of them may be rendered visible(quasi-invisibility). Finally in Section VI we summarize our results 
and conclusions. 


II. THE HIGGS SECTOR OF MINIMAL SUPERSYMMETRIC STANDARD MODEL 


The ATLAS [I, 2] and CMS [3, 4] experiments at the CERN LHC have independently observed a resonance, 
whose properties are consistent with the SM Higgs boson. The ATLAS experiment after collecting data at 
an integrated luminosity of 4.8 fb“^ at ^/s = 7 TeV in 2011 and 5.8 fb“^ at ^/s = 8 TeV in 2012 confirmed 
the evidence for the production of a neutral boson with a measured mass of 126.0 ± 0.4(stat) ± 0.4(sys) GeV, 
with a significance of 5.9 a. The CMS experiment after collecting 5.1 fb“^ at 7 TeV and 5.3 fb“^ at 8 TeV 
reported an evidence of a neutral boson at 125.3 ± 0.4(stat.) ± 0.5(syst.) GeV with a significance of 5.8 a. As 
mentioned in the Introduction, within the framework of the MSSM, we shall identify this resonance with the 
lightest Higgs boson of the model. We recall that the tree-level Higgs boson masses are determined by CP-odd 
Higgs boson mass Ma and tan/I. Requiring the production cross section of 126 GeV Higgs boson decaying to 
two photons agrees with the one observed at the CERN LHC, divides the MSSM Higgs parameter space into 
two distinct regions [26]: 

• The non-decouplng regime where Ma ^ 130 GeV. In this region the heavier CP-even state 77° is SM- 
like, and the light CP-even Higgs state /i° and the CP-odd Higgs state A° are almost degenerate in 
mass [27] and close to Mz, while the charged Higgs bosons are nearly degenerate with 77° [28]. The LHC 
phenomenology for this sector has been discussed in the past [29-31]. 

• The decoupling limit where Ma ^ 300 GeV. In this region the light CP-even Higgs boson /i° is SM-like, 
and all the other physical Higgs boson(s) are heavy and almost degenerate with A° [32]. The decoupling 
properties of the MSSM Higgs sector are not special to supersymmetry. Rather, they are a generic feature 
of non-minimal Higgs sectors. 

The non-decoupling scenario, which leads to a light SM-like Higgs, is highly constrained. Thus, the decoupling 
regime is a more viable scenario as far as the MSSM Higgs search results are concerned. Therefore, we shall 
consider only this scenario in this paper. For the non-decoupling regime of supersymmetric models, see, e.g. [25]. 

Having summarized the experimental results and the different scenarios, we now summarize the aspects of 
the MSSM Higgs sector that are relevant to our discussion. The masses Mi^/u and the mixing angle a of the 
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neutral CP-even Higgs states are well known. These can be written as [33] 

^h/H ~ 2 ^-^11 + ^-^22 T '\J— 2M^M'^C4i3 + C ), 

„ 2AM22 - [Ml + M|)s^ 

tj 9.11 Q; — -, , 

AM2^ - AMI 2 -b (M| - Ml)c2ii + 

where 


(11.1) 

( 11 . 2 ) 


C = 4AM]^2 + (AMfi - AM| 2 )^ - 2{Ml - M|)(AMi\ - AM|2)c2/3 - 4(Ml -b M|)AMf2S2/3, (11.3) 


and C 4 ^ = cos 4 / 3 ,C 2/3 = cos2/?,C /3 = cos/3, and S /3 = sin/3. Furthermore, AMij,i,j = 1,2 quantify the radiative 
corrections to the CP-even Higgs boson mass matrix. 

The dominant radiative corrections arising from the top-stop sector are contained in AM 22 . Identifying the 
resonance discovered at 126 GeV with /i°, one can, then, write the dominant radiative correction to the Higgs 
boson mass matrix in terms of the mass Mh of as 


2 M^{Ml 
AM^ = ^ 


Ml - Ml) - MlMlclp 




MXsl-Ml 


(H.4) 


where we have used the approximation AM 22 ^ AMn, AM 12 , and where Ma is the mass of the pseudoscalar 
Higgs boson We note that there is another solution for AMI 2 , which is unphysical [34], and is, therefore, 
not relevant. 

With this approximation we can write the mass of and the mixing angle a in terms of Ma and tan /3 as 


a 


(Ml + Ml- MDiMlcj + Mlsj) - MlMlcl^ 


— tan 


Mlcj + MW. - Ml 


{Ml- 


■ MDcpsp 


Mlcj + Mlsj-Ml 


(H.5) 

(H. 6 ) 


Thus, in principle one can calculate the mass of . However, in actual practice the mass of will depend on 
other parameters of the model, which include the supersymmetry conserving Higgs bilinear parameter /r, and 
the supersymmetry breaking scale Ms- In order to obtain a handle on the behavior of solutions of interest to 
us, we perform a numerical scan using the package CalcHEP [35]. For our study we use a set of inputs which are 
consistent with the known experimental constraints, and also which have the possibility of leading to spectra 
that may be visible in the near future. Thus, we are guided by the principle of search for SUSY in upcoming 
experiments and choices of parameters that continue to make low energy SUSY a viable option to address the 
naturalness and hierarchy problems of the SM. We have calculated the dependence of Mh on fj, for different 
values of tan /3 and the supersymmetry breaking scale Ms- In Figs. 1 and 2 we show the dependence of Mh on 
fjL for values of Ma = 300 GeV and Ma = 500 GeV. From these two figures we conclude that Mh does not vary 
significantly as a function of /x. For fixed values of tan /3 it has a weak dependence on M 5 . However, Mh has a 
significant dependence on Ma and tan /3, and can be described fully in terms of these two parameters when we 
use the fact that Mh is in the range 123-129 GeV. This feature had been pointed out earlier in Refs. [36-38]. 
Our results agree with these authors. In the present work these results have been obtained numerically by 
using as input the parameters Ms, tan/3 and ^ as detailed above, and adjusting the scalar trilinear coupling 
At so as to have Mh in the range 123-129 GeV. We have checked that the results of our calculations are in 
good agreement with those of [33]. 


III. THE NEUTRALINO SECTOR OF THE MSSM WITH GUT BOUNDARY CONDITIONS 

In order to study the invisible decays of to neutralinos, we consider the neutralino sector of the MSSM 
in some detail. For this purpose we consider the neutralino mass matrix of the model and the implications of 
the GUT boundary conditions on the neutralino spectrum. This analysis will lead to a general understanding 
of the nature of mixing between the gaugino and Higgsino states, as well as those of the physical neutralino 
states after electroweak symmetry breaking. 

We recall that the neutralinos are an admixture of the fermionic partners of the two Higgs doublets. Hi and 
H 2 , and the fermionic partners of the neutral gauge bosons. When the electroweak symmetry is broken, the 
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FIG. 1: Mh as a function of /x for Ma = 300 FIG. 2: Mh as a function of /x for Ma = 500 
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physical mass eigenstates are obtained from the diagonalization of the neutralino mass matrix [39-41] 


.M’mssm = 


Ml 

0 

—mz sin 9w cos f3 


0 

M 2 


mz cos 9w cos /3 


—mz sin 9 iy cos l3 sin sin/3\ 

mz cos 9w cos j3 —mz cos 9w sin /? 1 
0 -u 




(III.l) 


where Mi and M 2 are the G(l)y and the SU{2)l soft supersymmetry breaking gaugino mass parameters, /x 
is the Higgs(ino) mass parameter, Mz is the Z boson mass, 9\y is the weak mixing angle and tan/3 = U 2 /U 1 
is the ratio of the vacuum expectation values of the neutral components of the two Higgs doublet fields Hi 
and H 2 . The gaugino mass parameters Mi and M 2 may have some relation between them, or they can be 
completely arbitrary. If we assume that the MSSM is embedded in a grand unified theory, then the boundary 
conditions at the GUT scale imply a definite relation between the gaugino masses, which would imply a relation 
between them at the weak scale following the renormalization group evolution. Here we will consider two types 
of boundary conditions on Mi that follow from embedding of MSSM in a grand unified theory, namely the 
universal boundary conditions and the nonuniversal boundary conditions. 


A. Universal boundary condition 

In the MSSM, with universal gaugino masses at the grand unified scale, usually referred to as mSUGRA, the 
soft supersymmetry breaking gaugino mass parameters Mi and the corresponding gauge couplings ai{i = 1,2,3) 
satisfy the boundary condition 

Ml = M2 — Mu = mij2-, 

ai = a2 = q;3 = ac, (HI.3) 

at the grand unified scale Mg- Using the one-loop renormalization group equations [42] for the gaugino masses 
and the gauge couplings leads to the ratio 


Ml : M 2 : M 3 ~ 1 : 2 : 7.1, (HI.4) 

for the soft gaugino masses at the electroweak scale Mz- In the following, for definiteness, we shall consider the 
value of tan /3 = 10. If one is interested in Higgs boson decaying invisibly into light neutralinos, it is appropriate 
to consider a light or massless eigenstate of the mass matrix (III.l). It has been shown earlier [18, 25] that 
with the gaugino mass parameters satisfying the universal boundary conditions at the GUT scale, it is not 
possible to satisfy the masslessness condition for the lightest neutralino following from the determinant of the 
mass matrix (III.l) and the experimental constraints as implied by the LEP experiments [43]. Thus, a massless 
neutralino is ruled out with universal boundary conditions (HI.3) at the GUT scale. 
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TABLE I: Ratios of the gaugino masses at the GUT scale in the normalization Mi{GUT) = 1, and at the 
electroweak scale in the normalization Mi{EW) = 1 for F-terms in different representations of SU{5). These 
results are obtained by using 1 -loop renormalization group equations. 


B. Non-universal boundary condition 

Universal soft supersymmetry breaking gaugino masses are not the only possibility in a grand unified theory. 
In fact, non universal boundary conditions for the soft gaugino masses can arise quite naturally in SU (5), S'O(IO) 
and Eq supersymmetric grand unified theories. 

In grand unified supersymmetric models, gaugino masses are generated by a non-singlet chiral superfield 
that appears linearly in the gauge kinetic function /($), which is an analytic function of the chiral superfields 
<!> in the theory [44]. The gaugino masses are generated from the coupling of the field strength superfield IU“ 
with /(4>), when the auxiliary part F$ of a chiral superfield $ in /(4>) gets a VEV. When F$ gets a VEV 
< Fj> >, we obtain the Lagrangian containing the gaugino masses 

Cg.k. D + h.c., (III.5) 

IVl p 

where A®’** are gaugino fields. Here, we denote by A^, A^ and A^ the U{1)y, SU{2)l and SU{3)c gaugino 
fields, respectively. Since the gauginos belong to the adjoint representation of the gauge group, 4> and F$ can 
belong to any of the representations appearing in the symmetric product of the two adjoint representations of 
the unified gauge group. 

Eor example, in the case where the SM gauge group is embedded in the grand unified gauge group SU{5) 
the symmetric product of the two adjoint (24 dimensional) representations of SU{5) leads to 

(24 (g) = 1 0 24 0 75 0 200. (III. 6 ) 

In Table I we have summarized the ratios of gaugino masses which result when 0$ belongs to different repre¬ 
sentations of SU{5) in the decomposition (III. 6 ). 

Similarly, nonuniversal gaugino masses can arise from the embedding of MSSM in a grand unified theory 
based on S'O(IO) and Eq. For these gauge groups we have the decomposition 
50(10) : 


Ee : 


(45 0 45)symm = 1 0 54 0 210 0 770; 


(III.7) 


(78 0 78)symm = 1 0 650 0 2430. (III. 8 ) 

The analogs of Table I for the gauge groups SO(10) and Eq are given in Appendix A of [18], where a detailed 
discussion has been presented on the phenomenological consequences of different choices of the grand unified 
gauge group. We note that in the present work we are not necessarily looking at very light neutralinos as the 
second lightest Higgs boson can be heavy, and can decay into neutralinos that are massive. 


C. Composition of Neutralinos 


The composition of the lightest neutralino Xi in terms of the gauginos and Higgsinos can be written as [45, 46] 

X° = ZnB + Zi2fU3 0 013^0 0 014^0 (111.9) 
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FIG. 3: The contours of constant lightest 
neutralino mass m^o in MSSM with arbitrary 
gaugino masses at the GUT scale in the fj, — Mi 
plane. Here the value of the parameter M 2 = 200 
GeV. 



Mi(GeV) 

FIG. 4: The contours of constant second lightest 
neutralino mass m^o in MSSM with arbitrary 
gaugino masses at the GUT scale in the fj, — Mi 
plane. Here the value of M 2 = 200 GeV. 


where 


cot (Ml — m^o ) (Ml — m^o ) (/i sin (3 + m^o cos /3) 

M 2 — m^o ’ mz sin 9wim^o + sin 2/3) ’ 

The composition of the second lightest neutralino can be written as 



(Ml — m^o ) (/I cos 13 + m^o sin (3) 
mz sin dw (m^o + /i sin 2/3) 


(HI.IO) 


= Z 21 B + Z22W^ + ^23-^1° + ^24^2° 


(HI.ll) 


where 

(Ml — m^o) (/I cos /3 + m^o sin (3) 
mz sin 9w (m^o + /i sin 2/3) 

(HI.12) 

The possibility of decaying into neutralinos depends on the mass of the neutralinos into which it can decay. 
Our objective is to study the invisible decay of the second lightest Higgs boson into neutralinos. For this 
purpose, we carry out a study of the neutralino mass as a function of Mi and /i, as these parameters, apart 
from M 2 , control the neutralino mass matrix, and hence the neutralino mass eigenvalues. We have considered 
the constant mass curves for both the lightest and the second lightest neutralinos. Furthermore, we have 
considered each of the scenarios, namely the case of arbitrary gaugino masses, universal, and non-universal 
gaugino masses, respectively at the grand unified scale in Figs. 3, 4, 5, 6, 7 and 8. 

For arbitrary supersymmetry breaking gaugino masses at the grand unified scale, we shall choose the pa¬ 
rameter M 2 = 200 GeV, a choice we shall discuss in the following. It is obvious that choosing universal or 
nonuniversal boundary conditions at the grand unified scale is a rather restrictive condition. 

In Figs. 3 and 4 we plot the contours of constant lightest and the second lightest neutralino masses in the 
/r — Ml plane for arbitrary gaugino masses at the GUT scale. It is clear that m^o increases with Mi but does 
not change significantly with /x. On the other hand the second lightest neutralino mass m^o does not change 
appreciably with Mi but increases with /i. This is primarily because Xi is bino-like, and X 2 i® & Higgsino- 
like state. For low values of Mi, the mass of Xi is almost equal to Mi, and it is bino-like for all values of /i. 
Furthermore, when Mi > /i, x? is Higgsino-like. At large /x, and when Mi > M2(200GeV), Xi will be wino-like. 
We see from Fig. 3 that with arbitrary gaugino masses and with Mi larger than 150 GeV and /x larger than 
250 GeV the mass of Xi approaches 150 GeV and, thus, beyond these values a 300 GeV Higgs cannot decay 
into two lightest neutralinos. 

For low values of Mi and with /x < M2(200GeV), X 2 i® Higgsino-like. Furthermore, for values of /x > 
M2(200GeV), X 2 i® wino-like for all values of Mi. For values of /x < Mi, X 2 i® Higgsino-like. In Fig. 4 we see 
that for Ml > 150 GeV and > 200 GeV, the mass of X 2 approaches 150 GeV, and hence for values larger 
than these a 300 GeV second lightest Higgs cannot decay into two second lightest neutralinos. Thus, only when 



^ 2 ^ = 1 , - 


cot 9w (Ml — m^o ) (Ml — m^o) (/x sin (3 + m^o cos /3) 


X2^ 




M2 — m^o 
X.2 


mz sin 9w (?^xE 
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FIG. 5: The contours of constant lightest 
neutralino mass m^o in the /i — Mi plane in 
MSSM with universal soft gaugino masses at the 
GUT scale. 


FIG. 6 : The contours of constant second lightest 
neutralino mass m^o in the /i — Mi plane in 
MSSM with universal soft gaugino masses at the 
GUT scale. 



FIG. 7: The contours of constant lightest 
neutralino mass m^o in the fx — Mi plane with 
MSSM embedded in SO(IO), and with 
non-universal soft gaugino masses at the GUT 
scale. 



FIG. 8 : The contours of constant second lightest 
neutralino mass m^o in the fx — Mi plane with 
MSSM embedded in Eq, and with non-universal 
soft gaugino masses at the GUT scale. 


Ml is large and /i is relatively small, or if fx is large and Mi is relatively small, can decay invisibly into 
neutralinos. 

It is instructive to consider here the case of universal boundary conditions on the gaugino masses at the 
grand unified scale. Gonsidering the values of Mi in the range 200 — 400GeV, which implies values ofM 2 in the 
range of 400 — SOOGeV, for low values of fx, is Higgsino-like, and for large values of fx, where /r > Mi , fx, m^o 
is bino-like. Similarly, for values of M 2 > Mi,/x m^o is Higgsino-like. In Fig. 5 and Fig. 6 we plot the contours 
of constant neutralino masses in the /r — Mi plane for the case of universal boundary conditions on the soft 
gaugino masses at the grand unified scale. Here and m^o do not change significantly as a function of Mi. 
However, both the masses are an increasing function of fx. From Fig. 5 we see that with /i larger than 175 GeV, 
m^o approaches 150 GeV, and, thus, a 300 GeV Higgs cannot decay into two lightest neutralinos for the case 
of universal boundary conditions. For values of /i larger than 275 GeV m^o approaches 250 GeV for Mi > 300 
GeV. Hence, for these values of fx, a 500 GeV Higgs cannot decay into lightest neutralinos.. However, for low 
Ml larger values of fx will be allowed. In Fig. 6 we can see dominantly depends on fx. For fx larger than 
150 GeV, m^o crosses the 150 GeV limit and the decay of the 300 GeV second lightest Higgs into two second 
lightest neutralinos will be forbidden for these values of /i for any value of Mi. Similarly, for a 500 GeV Higgs 
boson, ^ = 250 GeV is the limit because for higher values of fx, m^o approaches a value of 250 GeV. 

In Fig. 7 and 8 we discuss two typical examples of nonuniversal boundary conditions on the gaugino masses 
at the GUT scale. The first is the case of SO(10) GUT with the condition that the three gaugino masses are 
in the ratio Mi : M 2 : M^=l : 6 : —14.3, and the second one is with Eq GUT with the boundary condition 
in which the three gaugino masses are in the ratio Mi : M 2 : M 3 = 1 : 50.2 : 70.9. The case of S'O(IO) is 
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akin to that universal boundary condtion discussed above because the range of M 2 values considered is greater 
than the range of Mi and The range of values of Mi considered here is 100 — SOOGeV. For Mi < fi, Xi is 
bino dominated, and for Mi > fi, Xi is Higgsino dominated. On the other hand, the particular choice of Eq 
representation is very interesting because this is the only GUT representation in which the lightest GP-even 
Higgs /i° can decay into two lightest neutralinos. With Eq boundary condition we can see that the mass of the 
lightest neutralinos does not depend on /i, whereas in the case of SO(IO) the mass depends on both Mi and fx. 
For this Eq representation the range of Mi values considered here is 20 — 80GeV, and values of Mi < M 2 
for the entire parameter space. Hence in this case Xi is always bino-like. 


We recall here that there are additional constraints coming from the LEP analysis of the decay in 
invisible modes, on the /i — Mi parameter space. The invisible decay of the lightest Higgs boson to the lightest 
neutralinos, if kinematically allowed, is mainly constrained by the Z invisible decay rate. This invisible decay 
width has been measured very precisely by the LEP experiments [43] with 

r{Z° ^ XiXi) < 3 MeV. (HI.13) 

The Z width to a pair of lightest neutralinos can be written as [47] 

( A 2 \ 3/2 

1-^) . (HI.14) 

m^o J 

For our analysis we have taken tan/3 = 10 [28]. The trilinear soft supersymmetry breaking coupling At has 
been adjusted in order to obtain a lightest Higgs boson mass of « 126 GeV. The gluino mass is taken to be 

1400 GeV [48], and the squarks are assumed to have a mass above 1 TeV [49], thereby respecting the current 

experimental bounds. 

In our calculations we have imposed the constraint of the lightest chargino mass bound of m^+ > 94 GeV 
following from the LEP experiments as well as the bound from invisible Z^ decay width coming from Z^ decay 
into neutralinos [50]. 

Having discussed in some detail the correlations among the parameters of the neutralino sector, we now turn 
to the implications of this analysis for the possible invisible decays of for a choice of parameters which are 
consistent with the constraints discussed in this section. 


IV. DECAY OF THE SECOND LIGHTEST HIGGS BOSON TO NEUTRALINOS IN THE 
MINIMAL SUPERSYMMETRIC STANDARD MODEL. 

We now address the main issue of this paper, namely the invisible decays of the heavier GP-even Higgs 
boson into neutralinos in the MSSM [51, 52]. We consider these decays in the decoupling regime of the 
MSSM, where the Higgs boson is relatively heavy. In the decoupling regime, the heavy can decay into 
neutralinos which are not necessarily light. 

For this purpose we catalog the decay widths of the heavier GP-even Higgs boson into a pair of lighter 
neutralinos in the MSSM. These can be written as 


r(M° x'j’xi) = - 4m|o/M|f)3/^ [(Z 12 - tan6»H/Zii)(Zi3COSQ;- Z^sina)]^ , (IV.l) 

2v27r 

where Zij are the elements of the matrix Z which diagonalizes the neutralino mass matrix (111.1). In the 
decoupling limit, when the mass of the pseudoscalar is large compared to the mass of the Z boson Mz^ the 
Higgs mixing angle a ^ f3 — 7 r/ 2 , so that the decay width (IV.l) can be written as 


r(if° ^ X?X?) = (1 - 4 m|o/M|,) 3/2 [(Z 12 - tan 0 u^Zn)(Zi 3 sin/3 + Z 14 cos/3)]' . (IV.2) 

2v27r 

Similarly, the decay width of EI^ into a pair of second lightest neutralinos can be written as 

A/f ^ A/f 

r(77° X 2 X 2 ) = ^ 2 vi 7 r ~ 4TO|o/M|f)3/^ [{Z 22 - tan 6 »h/Z 2 i)(Z 23 cos a - Z 24 sina)]' , (IV.3) 
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which in the decoupling regime becomes 

r(7J° ^ X 2 X 2 ) = (1 - 4w|o/M|f)^/^ [(Z 22 - tan6»M/Z2i)(Z23sin/3 + ^ 24003 /?)]^ . (IV.4) 


2v^71 


On the other hand the decay width for the process —>■ X 1 X 2 written as 


V27i 


X'i' 


x[l - - m\oJMjj - 2(eie2/M|r)m^oTO^o], 


where 


^"121 = 2 (^22 — tan0H^Z2i)(2'i3 cos a — Z14 sin a) + 2 (^12 — tan 04yZii)(Z23 cos a — Z24 sin a), (IV. 6 ) 

which in the decoupling regime reduces to 

_F'i2i = 2 (-^22 ~ t&a 9\y Z2i){Zi^ sin (3 + ^14 cos j 3 ) + —(^12 ~ t&a ^ 44)(^23 sin (3 + 2^24 cos /?), (IV.7) 


and where the constants CiS (i=l,2,3,4) stand for the sign of the neutralino mass. When the neutralino mass 
matrix is diagonalized, we allow the sign of the ith eigenvalue to be either positive or negative. 

In our previous work [18], we had considered the possibility of the lightest Higgs boson decaying into lightest 
neutralinos in the MSSM and the NMSSM. As the mass of the lightest Higgs boson is fixed to be 126 GeV, 
we had considered the case of the lightest neutralino with mass rn^o < mfiO/ 2 - We had concluded that such a 
constraint on the neutralino mass is not satisfied with either universal or nonuniversal boundary conditions on 
the soft gaugino masses at the GUT scale. The only exception to this conclusion was a nonuniversal scenario 
with a very large dimensional representation(2430) of the gauge group Eq. However, this is not an appealing 
possibility. The constraint on the light neutralino masses was obeyed in the MSSM only with arbitrary gaugino 
masses at the GUT scale. Since in the present paper we are considering the invisible decay of the heavier Higgs 
boson the lightest neutralinos need not be massless or very light. Thus, in the following we shall consider 
the case of massive neutralinos, which can arise in all the three cases, i.e. universal, nonuniversal and arbitrary 
soft gaugino masses at the GUT scale. 


V. RESULTS FOR THE INVISIBLE DECAY OF THE HEAVIER HIGGS BOSON 


Having summarized the results for the decay widths for the invisible decay of the heavier Higgs boson in the 
previous section, we now evaluate these decay widths using the parameter space of the MSSM allowed by the 
present experimental constraints. For this purpose, we shall use the boundary conditions as implied by the 
embedding of MSSM in a grand unified theory, as well as for arbitrary soft gaugino masses at the GUT scale. 
As discussed in the last section, the final states in the decay process that are of interest to us are 2. 

We first summarize the parameter space used in our analysis. The trilinear soft supersymmetry breaking 
parameter At pertaining to the stop is adjusted to obtain the lightest Higgs boson mass in the range 123-127 
GeV, although it does not affect the neutralino sector. To get the lightest Higgs mass in this range we use the 
’’maximal mixing” scenario [53], wherein 

At — pL cot j 3 = x/QMs, (V.l) 

where Ms is the soft supersymmetry-breaking scale, which we take to be 1.5 TeV. The gluino mass is taken 
to be 1.4 TeV, and squarks are assumed to have masses above 1 TeV. The Slepton masses are assumed to 
be greater than 500 GeV, thereby respecting all current experimental bounds [54]. For the case of arbitrary 
boundary condition on the gaugino mass parameters, we have taken M 2 = 200 GeV. Furthermore, we have also 
imposed the constraint of the lightest chargino mass ra^+ > 94 GeV as implied by the LEP experiments, as 
well as the bound on the invisible width into neutralinos. We note that decreasing the value of M 2 increases 
the chargino mass bound on /r, thereby eliminating the parameter space with large invisible branching ratio 
for the lightest Higgs boson [50]. 

In addition, we have also imposed the constraints resulting from (g-2) of the muon, as well as flavor constraints 
following from 6 —)■ 57 and Bg —> We have implemented these constraints using CalcHEP. The input 

parameters so obtained are summarized in Table H [55, 56] 
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tan/3 = 10 

Ms = 1.5 TeV 

Ma = 300,500 GeV 

Ma = 200 GeV 

Ms = 1402 GeV 

At= 3600 GeV 

Ai= 3600 GeV 

A.r= 1000 GeV 


TABLE II: Input parameters for the MSSM consistent with all the experimental constraints. 



Mi(GeV) 

FIG. 9: The contours of constant branching ratio 
for —>• XiXi) in the /r — Mi plane for 

Ma = 300 GeV with arbitrary gaugino masses at 
the GUT scale. Here M 2 is taken to be 200 GeV. 
The shaded region represents the region allowed 
by the WMAP data. 



FIG. 10: The contours of constant branching 
ratio for —>• XiAi) in the ^ — Mi plane for 

Ma = 500 GeV with arbitrary gaugino masses at 
the GUT scale. Here M 2 is taken to be 200 GeV. 
The shaded region is allowed by the WMAP 
data. 


We now proceed to present the results of our calculations. The results are presented as constant branching 
ratio contours in the p. — Mi plane for two different values of Ma ~ Mh = 300 GeV, 500 GeV. For the case of 
arbitrary soft supersymmetry breaking gaugino mass parameters we have taken M 2 = 200 GeV, and tan /3 = 10. 
In Fig. 9 we observe that the branching ratio for —)■ XiXi can be at most 16% for Ma = 300 GeV, and that 
too in a very narrow region of the parameter space. Furthermore, from Fig. 10 we can see that this branching 
fraction reduces with the increase in the value of Ma, having a value of 10% for Ma = 500GeV. 

Since we are considering MSSM with Rp conservation, the lightest neutralino is an absolutely stable particle. 
In this case it is important to check that it is not overproduced in the early universe. We have, therefore, 
calculated the relic density of the lightest neutralino for the range of parameters used in our calculations using 
micrOMEGAs [57], and imposed the constraints following from the WMAP data [58, 59]. The result of this 
calculation is shown in Fig.9 as the shaded region allowed by the constraints of the relic density considerations. 
Similarly, in Fig. 10 the shaded region is the region allowed by the WMAP data. 

This behavior of the decay —)■ XiXi iii&y be understood as follows. At low values of Mi, x? is bino-like 

for all values of /i in the region of the parameter space that we are considering. The Higgs decay to neutralino 
is suppressed because of small value of Zn for tan /3 = 10. The mass of the lightest neutralino increases with 
the increase in the value of Mi. In Fig. 9, for Mi > 150 GeV there is a kinematic suppression for Ma = 300 
GeV, resulting in a smaller branching ratio. On the other hand for Ma = 500 GeV there is no kinematic 
suppression, hence the branching ratio does not fall off, as is clear from Fig. 10. 

For small values of Mi and for all values of fi in the parameter space that we are considering, the quantity 
(Ml — TO^o) is close to zero, so that the decay width for H —> XiXi is reduced as can be seen from the analytical 
results in the previous Section, Eq. (IV.3). From Eq. HI. 10, we observe that Z 13 is almost zero. Hence the 
branching fraction for both values of Ma = 300 GeV and 500 GeV is small. For Hxed /i, with increasing 
Ml > fi, Xi is Higgsino-like, (Mi — m^o) increases, the total width decreases, and hence the branching ratio 
increases. For fixed Mi with decreasing /i, (Mi — m^o) increases. Furthermore, the denominator (total width) 
decreases with increasing /i, the numerator decreases faster, and hence the branching fraction decreases with /i. 
Hence, to obtain a constant branching fraction, we have to increase Mi and /i simultaneously until the factor 
(1 — 4m|o/M^)^/^ starts dominating, and causes kinematic suppression for the case of Ma = 300 GeV. 
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Mi(GeV) 

FIG. 11: The contours of constant branching 
ratio of {H^ —> X 1 X 2 ) the ^ — Mi plane for 
Ma = 300 GeV with arbitrary gaugino masses at 
the GUT scale. Here M 2 is taken to be 200 GeV. 
The shaded region is allowed by WMAP data. 



Mi(GeV) 

FIG. 12: The contours of constant branching 
ratio for {H^ —)• X 1 X 2 ) the — Mi plane for 
Ma = 500 GeV with arbitrary gaugino masses at 
the GUT scale. Here M 2 is taken to be 200 GeV. 
The shaded region is allowed by WMAP data. 


We note that the total width has been calculated using CalcHEP, and depending on the values of parameters 
Ml and fi different channels contribute to the total width, which has been used in the calculations presented 
here. 

We may mention here that if we increase the value of M 2 > 200GeV, large branching fractions for the 
invisible decay would still be possible. On the other hand, values of M 2 < 200GeV will exclude the region of 
parameter space with considerable invisible branching fraction. 

We now turn to the invisible decay —)■ X 1 X 2 values of Ma ~ Mh = 300 GeV and 500 GeV. In Fig. II 

and Fig. 12 we show the contours of constant branching ratios for this decay. For Ma = 300 GeV the branching 
ratio can at most be 12% and for Ma = 500 GeV it is at most be 6%. As we have seen in case of the decay 
XiXi s-lso we see a kinematic suppression in case of 300 GeV Higgs boson for large values of Mi 
and /X. The reason is that the sum of the masses of the lightest and the second lightest neutralino reaches its 
limiting value and the decay —> X 1 X 2 longer possible. On the contrary, for 500 GeV Higgs boson there 
is no kinematic suppression, because M^ which is nearly equal to Ma in the decoupling regime is sufficient to 
produce a lightest and a second lightest neutralino for this particular region of the parameter space and the 
decay width is governed by the function F 121 in Eq.( IV.7). We note that the shaded region in these Figs, 
represents the region allowed by the WMAP data. 


In Fig. 13 and Fig. 14 we show the contours of constant branching ratio for the decay —> X 2 X 2 f^e 
/i — Ml plane for values of Ma ~ Mh = 300 and 500 GeV, respectively. The branching ratio can at most be 
2% in the case of 300 GeV second lightest Higgs, and it reaches a value of about 5% in case of 500 GeV second 
lightest Higgs boson. There is a kinematic suppression for 300 GeV Higgs for values of fi larger than 200 GeV as 
discussed earlier but for 500 GeV Higgs there is no such suppression and the branching ratio actually increases 
with increasing Ma in some regions of the parameter space. The shaded region in these Figs, represents the 
region allowed by the WMAP data. 

The mass of the second lightest neutralino becomes almost equal to Mi for large values of fx and Mi. 
Hence, in this region the branching ratio goes to zero. And from the Fig. 4 we can see that for values larger 
than Ml = 100 GeV, m^o does not depend on Mi and depends only on Hence, for a fixed values of fi, 
the quantities Z 22 , .^23 and Z 24 increase with Mi as (Mi — m^o) increases with Mi. Therefore, the width 
X 2 X 2 also increases. We have checked that the total width decreases with Mi. So the branching fraction 
> X 2 X 2 increases with Mi. Now for a fixed Mi, the quantity (Mi — m^o) decreases with /i. Furthermore, 
Z 22 decreases with /i. Z 23 and Z 24 also decrease, having (Mi — m^o) in the numerator and irn ^2 +yxsin2/?) in 
the denominator. Also, the total width decreases with /x, but the numerator in the branching ratio decreases 
much faster. Hence the branching fraction —)■ X 2 X 2 decreases with /x. Let us recall the branching ratio 

increases with Mi. Gonsequently we have to increase both Mi and /x to get the contours of constant branching 
ratio. This is reflected in Figs. 13 and 14. 
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FIG. 13: The contours of constant branching 
ratio for (iJ° —> X 2 X 2 ) in. the /i — Mi plane for 
Ma = 300 GeV with arbitrary gaugino masses at 
the GUT scale. Here M 2 is taken to be 200 GeV. 
The shaded region represents the region allowed 
by the WMAP data. 
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FIG. 15: The contours of constant branching 
ratio for (iJ° —> XiXi) in the /i — Mi plane for 
Ma = 500 GeV with universal gaugino masses at 
the GUT scale. The shaded region is allowed by 
WMAP data. 



FIG. 14: The contours of constant branching 
ratio for {H^ —)• X 2 X 2 ) the fx — Mi plane for 
Ma = 500 GeV with arbitrary gaugino masses in 
the — Ml plane. Here M 2 is taken to be 200 
GeV. The shaded region is allowed by WMAP 
data. 



Mi(Gev) 

FIG. 16: The contours of constant branching 
ratio for —)• X 1 X 2 ) /i — Mi for 

Ma = 500 GeV with universal gaugino masses at 
the GUT scale. The shaded region is allowed by 
WMAP data. 



Mi(GeV) 

FIG. 17: The contours of constant branching ratio for —)■ X 2 X 2 ) the fi — Mi plane for Ma = 500 GeV 

for universal gaugino masses at the GUT scale. The shaded region is allowed by WMAP data. 





13 




FIG. 18: The contours of constant branching 
ratio for (iJ° —> XiXi) in th® M ~ -^i plane for 
Ma = 500 GeV for SO(10) with non-universal 
gaugino mass boundary condition. The shaded 
region is allowed by WMAP data. 


FIG. 19: The contours of constant branching 
ratio of {H^ —)• XiXi) in the — Mi plane for 
Ma = 500 GeV for Eq with non-universal 
gaugino mass boundary condition. The shaded 
region is allowed by WMAP data. 


Having considered the case of arbitrary gaugino masses at the GUT scale in detail, we now turn to the case 
of universal boundary conditions. In Figs. 15, 16 and 17 we have plotted the contours of constant branching 
ratio for —)• XiXi j X 1 X 2 : —)■ X 2 X 2 J respectively for a 500 GeV second lightest Higgs boson 

with universal boundary condition on the gaugino masses at the grand unified scale. In this case the boundary 
conditions imply M2 « 2 Mi and M3 « 7 Mi at the electroweak scale. Here we can see the invisible branching 
ratios for XiXi can attain a value 9% and X 1 X 2 can attain a value of 10%. The branching ratio —>• X 2 X 2 can 
at most be 2%. We have taken the range of Mi to be 200 GeV to 400 GeV for this case. For this case M 2 and 
M 3 depend on Mi and to respect the experimental constraints we need to take M 3 i.e the gluino mass above 
1400 GeV. As discussed in the previous section for values of ^ larger than 275 GeV, the mass of the lightest 
neutralino exceeds the limiting value of 250 GeV. Hence, a 500 GeV Higgs boson cannot decay into lightest 
neutralino pairs. In the case of the second lightest neutralinos, /i = 250 GeV is the limiting value because the 
mass of the second lightest neutralino attains a value of 250 GeV as can be seen from Fig. 6 . 

We note from the behavior of the branching ratios that m^o and do not have a significant dependence 
on Ml but depend on /i. This pattern is depicted in Figs. 5 and 6 . The behavior of —)■ X 2 X 2 can also be 

understood in the same manner. We note that m^o changes linearly with fj, and in large parts of the parameter 
space it is practically equal to /i. Hence from Eq. HI.12 we can see that the width of —>■ X 2 X 2 increases with 
Ml. We have checked that the total width decreases with increasing Mi. Hence the branching fraction also 
increases with Mi. The total width decreases with /i. Hence to get constant branching ratio we need increasing 
values of Mi with decreasing values of /i. This can be seen from Fig. 17. The behavior of —> XiXi is also 
similar to that of —> X 2 X 2 - The branching fraction is dominated by (Mi — m^o) term which increases with 
Ml. We note that for large values of Mi the behavior is almost same as that of —)■ X 2 X 2 ' The shaded region 
in these Figs, represents the regions allowed by the WMAP data. As we can see from Figs. 15, 16 and 17 
the universal boundary condition scenario is very much constrained by the WMAP data. Only small regions 
of the parameter space are allowed which can give rise to some invisible branching ratio. 

In Fig. 18 we have considered an example of non-universal boundary condition on the gaugino masses. For 
this we have chosen a 210 dimensional representation of the gauge group SO(IO) in which Mi : M 2 : M 3 is in 
the ratio 1:6: —14.3. We have plotted the contours of constant branching ratio for 77° —>■ XiXi the /x — Mi 
plane. We can see that the branching ratio can go upto 10% for a small region of the parameter space. The 
shaded region in this Fig. represents the region allowed by the WMAP data. 

In Fig. 19 we have considered another example of non-universal boundary condition on gaugino masses for 
a 2430 dimensional representation of the gauge group Eq. The gaugino masses are in the ratio Mi : M2 : M 3 
= 1 : 50.2 : 70.9. This particular representation is somewhat special. Because out of all the non-universal 
boundary conditions this is the only representation in which r(defined as the ratio Mi/M2) satisfies the the 
condition r < 0.04, resulting in a light neutralino. This actually allows the 126 GeV Higgs to decay into lightest 
neutralinos. In this case we can see the branching ratio for 77° —>■ XiXi can go upto 11% for a tiny region of 
the /i — Ml parameter space. The shaded region in this Fig. represents the region allowed by the WMAP data. 

For the case of Fig. 18 we see that for a particular value of Mi the branching fraction decreases with /i. 
Because the quantity Mi —decreases with /i at a particular value of Mi, increases with /x as can be seen 
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tan tan j8 

FIG. 20: Branching ratio of —>■ XiXi) with FIG. 21: Branching ratio of —?► X 1 X 2 ) with 

Mh = 300 GeV(red),500 GeV(blue) as a = 300 GeV(red),500 GeV(blue) as a 

function of tan/?. function of tan/?. 



tan p 

FIG. 22: Branching ratio of {H° —)■ X 2 X 2 ) with Mh = 300 GeV(red),500 GeV(blue) as a function of tan^. 


from Fig. 7. Hence Z 13 and Z 14 decrease with /i. The term Z 12 also decreases with /i, because the denominator 
6 M 1 — TO^o decreases with /i, but the rate of decrease is lower than the rate of increase of Mi — rrij^o in the 
numerator. The total width decreases with /i, but here also the rate of decrease is less, hence the branching 
fraction decreases with ^ for a particular value of Mi. 

Fig. 19 corresponds to the situation where the mass of the lightest neutralino is almost equal to Mi for all 
fjL, as can be seen from Fig. 8 . Hence {Mi — m^o) is a very small quantity which is more or less constant as a 
function of fx. Hence ^ 12,^13 and Z 14 all decrease with Mi for a particular having m^o, which is equal to 
Ml, in the denominator. But the total width also decreases with Mi. Hence the branching fraction remains 
almost constant as a function of Mi. 

It is known that the 126 GeV Higgs boson of the MSSM can decay invisibly only for arbitrary gaugino 
masses at the GUT scale. Here we find that for the situation for the second lightest Higgs boson MSSM is far 
more interesting, because it can decay invisibly, not only for arbitrary gaugino masses but also for constrained 
boundary conditions of universal and non-universal gaugino masses at the GUT scale. 

From the plots we notice that there are kinks in all the contour plots for branching ratios studied in this 
paper. Although it is difficult to obtain analytical expressions for the constant branching ratio contours, one 
thing that can be easily understood is that depending on the mass of the neutralinos different decay channels 
will open up for different range of values of /r and Mi. Whenever the total decay width in the denominator 
changes because of new channels opening up, the constant branching ratio contour plots show kinks. 

The range of values of Mi and considered in this paper is adequate for our study, because beyond this 
range of Mi and /i the branching ratios fall and are therefore not relevant. 


We have also considered the dependence of the branching ratios for the decays 77° 


XiXi 




X 1 X 2 


and 77° —>■ X 2 X 2 tan/3. Figs. 20, 21 and 22 show the dependence of these two branching ratios on tan/3 for 
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M^CGeV) 


FIG. 23: The branching fraction of the heavier Higgs to bb(gieen), neutralinos(red), charginos(black), rr 

(blue), lightest Higgs(orange) as a function of Mff. 


Ma ~ Mh= 300 and 500 GeV for the choice of parameters shown in Table III. These values of parameters are 
consistent with all current experimental constraints. 


Ml = 120 GeV 

Ms = 1.5 TeV 

Ma = 300,500 GeV 

Ma = 200 GeV 

M 3 = 1402 GeV 

At= 3600 GeV 

Ab= 3600 GeV 

200 GeV 


TABLE III: Input parameters for the plots of branching ratios versus tan /3. 

The behavior of the tan /3 dependence of the invisible branching ratios can be understood in the following 
manner. In MSSM, the bb coupling is times its SM value, and —>• tt coupling is times the 

SM value. 

At low values of tan /3 the coupling to the up-type quarks is large, hence in this region its decay to 
top-antitop quarks is dominant, if kinematically allowed. If the second lightest Higgs has a mass of 500 GeV, 
the branching ratio for this decay is 57%. In this region, the decay channel bb or the invisible decay channels 
are not signihcant. If the Higgs is not so heavy (around 300 GeV), then top quark channel is kinematically not 
available, and the main decay channel then will be to two lightest Higgs and some of the invisible channels. 

In case of moderate tan /3, say tan /? = 7, invisible branching ratio for —)■ XiXi peaks for both 300 GeV 
as well as 500 Gev Higgs boson. 

The branching ratios for 77° —>■ x°X 2 peaks at moderate values of tan/3 for 300 and 500 GeV second 
lightest Higgs. On the other hand X 2 being very heavy, the 300 GeV second lightest Higgs cannot decay into 
two X 2 S for any tan /3 for the parameter space considered. For the 500 GeV second lightest Higgs the branching 
ratio for 77° —>• X 2 X 2 peaks at moderate values of tan/3. 

For large values of tan ,3 « 30, Higgs coupling to the down-type quarks is dominant, and bb becomes the 
dominant decay mode. The invisible decay channels are once again insignificant. 

In the case of lightest Higgs decay, in the low tan (3 region, the invisible decay would be significant as the top 
decay channel is kinematically closed. In the case of the second lightest Higgs boson, low values of tan/3 do 
not give significant invisible decay width, but for moderate values of tan/3 « 10 the invisible branching ratio 
can be significant. 

In the Fig. 23 we have considered the branching fraction of heavier Higgs boson to all the dominant decay 
channels. The dominant visible channels are bb, tt, hh [60, 61]. With increase in the values of Mh, the branching 
fraction for the bb channel decreases, tt, hh channels are not significant, and chargino and neutralino channels 
become important. Neutralino branching fraction can be as large as 30%, and chargino branching fraction can 
be about 50 %. For this study the input parameters are Mi = 150, M 2 = 200, M 3 = 1400, = 200 and 

tan/3 = 10. From this study we can see that if the second lightest Higgs is heavy enough it can have major 
decay channels into the electroweak-inos. It can be inferred that in the future LHG run during the direct search 
for the second lightest Higgs boson, its decay into supersymmetric particles may play an important role in the 
decoupling regime. 
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A. “Quasi-invisible” decays 

We would like to draw attention of the reader to the fact that all the decay channels of the second lightest 
Higgs boson that we have considered and studied in detail in this paper are not truly invisible. —)■ XiXi 

channel is truly invisible as the lightest neutralino Xi is the lightest supersymmetric particle. It is the only 
stable supersymmetric particle assuming conservation of R-parity. All other heavy neutralinos will decay into 
the LSP and SM particles. Hence > X 1 X 2 > X 2 X 2 cannot be considered invisible. In ref [62] 

the authors have extensively studied the electrowikino decays for different regions of the parameter space of 
the electroweak sector without assuming any SUSY-breaking mediation scenario. There is a region of the 
parameter space defined as scenario C in the paper where j/rj < Mi, M 2 . In this region of the parameter space 
Xi and X 2 both Higgsino like and their masses are almost equal to the Higgsino mass /i and hence they 
are almost mass-degenerate. It has been pointed out in the same reference that the LSP multiplet production 
will be difficult to observe at the LHC because of the mass degeneracy and the soft decay products especially 
when the production is suppressed. Hence one has to be extremely careful while looking at the decay of the 
electrowikinos at the future LHC. Where the Xi and X 2 are almost mass degenerate, we have termed this kind 
of decays —)• X 1 X 2 —)• X 2 X 2 as ’quasi-invisible’ decay. It is worth mentioning that, one must expect 

that the International Linear Collider(ILC) will be able to identify the soft decay products, leptons and jets of 
mass 10 GeV or less, which are produced in the mass-degenerate case, because ILC has a cleaner environment 
for event reconstruction [62]. 

Another intriguing possibility that could render some of these channels “visible” is to adapt a strategy that 
has been used in the past for dark matter searches. In other words a promising search for invisible decays is in 
the monojet channel as in the case for dark matter search at the LHC. When the second lightest Higgs boson 
decays into invisible particles, monojet searches involving initial state radiation from the gluons that fuse to 
produce the Higgs boson could be used in these channels so the the Higgs decaying into light neutralinos does 
not necessarily escape the detection at the colliders. This same procedure is used to look for dark matter at 
the LHC. In case of invisible decays one has to look for monojet signals —)■ XiXi + jet as Xi is the lightest 
supersymmetric particle . In ref [63-65] it has been shown that the monojet signature carries a good potential 
to constrain the invisible decay width of a « 125GeV Higgs boson in a model independent fashion using the 
monojet search results by ATLAS and GMS [66, 67]. Now in case of MSSM with two GP-even Higgs, we can 
employ a different strategy while looking for the second lightest Higgs boson. One can calculate the production 
cross section x invisible branching in the MSSM for specific regions of the parameter space. There are several 
Standard Model processes which can act as background for the monojet signals, pp —)■ Z{-^ vv) + jet is the 
main irreducible background with the same topology as our signal. There is QCD background as well. The 
backgrounds can be estimated [68]. The background can be reduced significantly using several pr cuts on the jet 
and missing transverse energy cut. Then one can find the dependence of signal significance S=Ns/\/N b{N s is 
the number of signal events and Nb is the number of background events) on different parameters for 14 TeV HL- 
LHG with the desired integrated luminosity C. One can use the LHG monojet search results at 14 TeV, ie. the 
limits on the monojet events to probe the regions in the MSSM parameter space spanned by M^, tan/3 ,p [69]. 
Thus one can exclude regions of the MSSM parameter space at a desired signal significance(90 % or 95 % C.L). 


VI. DISCUSSION AND CONCLUSION 

In this work we have considered the possibility of invisible decays of the second lightest Higgs boson(M*') 
in the MSSM in the decoupling regime. In the past, various studies have shown that certain regions of the 
parameter space of MSSM, allow a Higgs boson in the mass range 123-129 GeV both in the decoupling and 
non-decoupling regime, satisfying the LHG constraints. For most of the parameter space, the lightest Higgs 
decay to the lightest neutralinos is kinematically allowed, leading to invisible decay modes. The main objective 
of those works was to prove that it would therefore be very important to study the couplings of the newly 
discovered particle at high precision. Global fits have been performed on the couplings of the newly discovered 
particle, in order to place upper bounds on the invisible decay width. Taking into account these bounds, the 
parameter space of these new physics scenarios can be further constrained, since the regions giving a large 
invisible Higgs decay branching ratio will be in conflict with the experiments. 

In the present work, which is a sequel to previous investigations [18], we are looking at the problem from a 
different point of view. We are asking a related question. We have considered the intriguing possibility that the 
heavier GP even Higgs(i7°) of the MSSM would decay into invisible LSP’s. Knowing the invisible branching 
ratios of the second lightest Higgs boson one can look for it in the future HL-LHG or ILG through the monojet 
signals. From the limits given by the colliders on the monojet signal one can then probe the MSSM parameter 
space at a desired confidence level. Hence it will also be possible to detect or rule out the second lightest 
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Higgs boson in the MSSM in a certain mass range from the monojet searches for its invisible decay. In other 
words, in our previous paper [18], the study has been done to constrain the MSSM parameter space looking 
at the possible invisible decay width of the newly discovered 126 GeV Higgs boson. But in this work we are 
following a different route and trying to make a theoretical prediction on the search for the second lightest Higgs 
boson of the MSSM. Our point is to understand the invisible decay width of consistent with the current 
experimental constraints. In the scenarios that we considered earlier, which partly motivated the present study, 
we had considered that the 125 GeV Higgs boson had some partial decay channels into light neutralinos in the 
MSSM, or that it would decay into neutralinos or light GP odd Higgs particles which could be present in the 
spectrum of the NMSSM. In [25] the intriguing possibility was considered that the 125 GeV Higgs was not that 
lightest GP even Higgs but the heavier one, while the lightest evaded detection altogether. Here, we have more 
conservative assumptions, viz., that the 125 GeV resonance is indeed the lightest CP even Higgs, but consider 
the possibility that the heavier one has invisible decays. This scenario can arise naturally in the decoupling 
region, i.e., the second lightest Higgs boson is quite heavy and could be produced at future experiments at 
the LHC which will eventually decay into lighter particles. These final states include neutralinos, charginos, 
bb, ZZ,TT,hh. In the present work we have tried to carry out an exhaustive search of the MSSM parameter 
space where the second lightest Higgs boson can decay into the lightest neutralinos and therefore remains 
invisible. We have considered the decay of the second lightest Higgs into the second lightest neutralinos also, 
because there are certain regions of the parameter space where the lightest and the second lightest neutralino 
are almost mass-degenerate and hence the second lightest neutralino may remain invisible(‘quasi-invisible’) at 
the LHC. We have used semi-analytical formulas as a guide to our study. In the decoupling regime, which is a 
realistic scenario, the values of Ma and Mh are nearly equal. 

For our parameter choices guided by recent constraints by LEP and LHC, we have taken the gaugino masses 
at the electroweak scale to be within present constraints. We have chosen At judiciously to get the lightest 
Higgs mass in the range in which it has been detected. We have taken soft SUSY breaking scale Ms to large, 
around 1.5 TeV consistently. In this paper, we have presented a detailed analysis first of the Higgs sector as well 
as that of the neutralino sector, in order to isolate the conditons under which it becomes kinematically feasible 
for (iL°) to decay invisibly. We have considered three different scenarios and have scanned the parameter 
space in terms of invisible decay of the second lightest Higgs boson in these three scenarios, namely arbitrary, 
universal and non-universal boundary condition on the soft gaugino masses at the GUT scale. Our analysis 
reveals that it is typically not possible to have such regions in the MSSM with GUT scale universal boundary 
conditions on the soft gaugino masses. This can be seen from Figs. 15, 16 and 17. The situation remains 
similar even with non-universal boundary conditions as can be seen from Figs. 18 and 19. The main reason for 
this is that there is not sufficient freedom in the choice of the gaugino masses with these boundary conditions. 
On the other hand, relaxing all constraints, as is the case with general (arbitrary) boundary conditions, allows 
the invisible branching ratios to be considerable. This can be seen from Figs. 9, 10, 11, 12, 13 and 14. 

From our study we conclude that there is a significant portion of the parameter space where the invisible 
decays can be quite significant. For instance, with universal or non-universal boundary conditions, the invisible 
branching ratio is not enhanced, but if we relax this constraint we can have more significant branching ratio in 
invisible decays. Hence the monojet searches for the invisible decay of the second lightest Higgs boson at the 
future colliders can be used as a probe to look for the second lightest Higgs boson as well as to put constraints 
on the regions of the MSSM parameter space. This would be most useful in the case of arbitrary boundary 
condition on the gaugino masses at the GUT scale. One has to be more careful and find different strategies in 
the case of universal and non-universal boundary condition. 
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